The A y Problem for p — 3 He Elastic Scattering 
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We present evidence that numerically accurate quantum calculations employing modern internu- 
cleon forces do not reproduce the proton analyzing power, A y , for p — 3 He elastic scattering at low 
energies. These calculations underpredict new measured analyzing powers by approximately 30% at 
Ec.m. = 1-20 MeV and by 40% at E c . m . = 1.69 MeV, an effect analogous to a well-known problem in 
p-d and n-d scattering. The calculations are performed using the complex Kohn variational principle 
and the (correlated) Hyperspherical Harmonics technique with full treatment of the Coulomb force. 
The inclusion of the three-nucleon interaction does not improve the agreement with the experimental 
data. 
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Over the past decade, numerical calculations for three- 
nucleon systems have reached a high degree of precision. 
This has made it possible to carry out accurate quan- 
tum mechanical computations for a variety of processes 
including 1) N-d elastic scattering and breakup 0J§], 2) 
N-d radiative capture 3) photodisintegration of 3 H 



and 3 He 
tering 



D , and 4) electron- 3 H and electron- 3 He scat- 
(see Ref. || for a more complete list of ref- 
erences). The calculations use a variety of approaches 
||, such as the Kohn variational principle, the Green 
Function Monte Carlo method, or direct solution of the 
Faddeev equations, and have made it possible to test our 
knowledge of the pairwise nucleon-nucleon (NN) interac- 
tion and to study effects of possible three-nucleon (3N) 
forces. 

Calculations employing modern NN and 3N interac- 
tions are capable of describing most of the experimental 
results for the processes listed above. However, there 
is a well-known and large discrepancy for the N-d an- 
alyzing power at low energies. What one finds is that 
the predicted A y values are sigificantly smaller in magni- 
tude than the measurements for both p-d and n-d elastic 
scattering. Resolving this "Ay puzzle" is a current and 
important area of research jl]|J . 

In the present Letter, we extend some of the analy- 
ses described above to the four-nucleon system. We will 
present new accurate computations for p- 3 Re elastic scat- 
tering at low energies. We will also report new measure- 
ments of the analyzing power for p- 3 He scattering which 
have been obtained for the purpose of testing whether 
the A y problem occurs in this system as well. 

Extending the accurate quantum calculations into the 
AN system is obviously of importance since it allows 
many new and stringent tests of the nuclear interaction. 
For example, accurate calculations of the alpha particle 
binding energy £? 4 have been achieved recently [^o| p^[ , 
and it has been found that the experimental value of the 



binding energy is reproduced with the same NN and 3N 
potentials that fit the 3 H binding energy. Since it appears 
that no four-nucleon potential is necessary to explain the 
ev-particle binding energy, one might expect that NN + 
3N interactions alone would also be sufficient to describe 
the various four nucleon scattering processes. 

The development of techniques for solving 4N prob- 
lems is also important for other reasons. Many reac- 
tions involving four nucleons, like d + d — > 4 Hc + 7 or 
p + 3 Hc — ► 4 He + e + + v e (the hep process), are of ex- 
treme interest from the astrophysical point of view. The 
theoretical description of these processes constitutes a 
challenging problem from the standpoint of nuclear few- 
body theory. Its difficulty can be appreciated by consid- 
ering for example the hep process. In Ref. fl3|| it was 
found that the capture from the initial P-wave channels 
("forbidden" transitions) gives about 40% of the calcu- 
lated S- factor, and this fraction depends critically on the 
correct description of the dynamics of the continuum and 
bound 4N states. 

Moreover, the study of 4N systems is important also 
for testing the various many-body techniques developed 
for studying systems with large (> 4) numbers of par- 
ticles. These theories necessarily include a number of 
approximations whose consequences can be investigated 
by comparison with more sophisticated calculations. The 
4N system is the simplest system in which these checks 
can be advantageously performed. 

The new calculations reported in this Letter are based 
on the Kohn variational principle and make use of cor- 
related hyperspherical harmonic (CHH) functions. The 
same approach has already been applied successfully for 
a variety of 3N processes [||JlJ]. A previous application 
for studying 4N scattering at zero energy was already re- 
ported [jlq] . In the present paper, we have improved our 
calculation to determine P- and D-wave phase shifts. 
The convergence of the P-wave shifts is rather slow and 
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has required a large technical effort to be achieved. Re- 
sults will be reported for calculations based on 1) the 
Argonne V18 NN potential Q (the AV18 model), and 
2) the Argonne V18 NN potential plus the Urbana IX 3N 
potential Q (the AV18UR model). The present calcu- 
lations represent the first attempt to study the effects of 
3N forces on the p- 3 He scattering observables at energies 
greater than zero. In this paper we focus on low energies 
where the convergence properties of our theoretical ap- 
proach are more satisfactory and where meaningful com- 
parisons with the experimental data can be performed. 

The new measurements of the proton analyzing power 
A y for p- 3 He elastic scattering were obtained at E c , m , = 
1.20 and 1.69 MeV. The experiments were carried out at 
the University of Wisconsin tandem accelerator labora- 
tory. Polarized protons from a crossed-beam polarized 
ion source fl7j| were accelerated, momentum analyzed by 
a 90° bending magnet, and transported to a 1-m scatter- 
ing chamber. The scattering chamber was filled with 43.4 
Torr of 99.95% purity 3 He gas, and was isolated from the 
beamline vacuum by a 4.44 x 10~ 5 cm thick Ni entrance 
foil located 1.27 cm from the chamber center. 

Elastically scattered protons were detected by three 
rectangular silicon surface-barrier detectors, 60 to 100 
/im thick, placed symmetrically on each side of the scat- 
tering chamber. A slit assembly restricted the angular 
field of view to ±0.34°. The spectra were clean except for 
a small contaminant peak that was well-separated from 
the peak of interest except at the most forward angle. At 
that angle, background subtraction was performed. 

After passing through the scattering chamber, the 
beam entered a polarimeter in which the beam polariza- 
tion was determined using p-a elastic scattering |l§| ] . The 
polarimeter was filled with one-half atmosphere natural 
He gas, and was separated from the scattering chamber 
by a 2.54 x 10~ 4 cm thick Havar foil. Because of the low 
beam energies, we could not measure the beam polariza- 
tion very accurately at the same time as data were being 
taken. However, previous experience indicates that the 
beam polarization does not normally change significantly 
over time. Consequently, at least once a day we made 
a careful measurement of the beam polarization by in- 
creasing the beam energy to 4.0 MeV at the center of the 
polarimeter. At this energy, the polarimeter analyzing 
powers are known to 2%. Each such careful measure- 
ment of the beam polarization yielded a value between 
0.79 and 0.84 with typical statistical errors of ±0.01. 

The new measurements are shown in Fig. [lj The er- 
ror bars include statistical uncertainties and also at the 
extreme forward angle an estimate of uncertainty in the 
background subtraction. There is also a scale factor un- 
certainty of 3% due to beam polarization uncertainties. 

We turn now to the theoretical calculations. Four- 
nucleon scattering problems have been studied theoret- 
ically for a long time (see Ref. |H| and references cited 
therein). Recently, increases in computing power have 



opened the possibility for accurate calculations of the 4N 
observables using realistic models for NN forces. These 
calculations have been performed mainly by means of 
the Faddeev-Yakubovsky (FY) approach |2^,|l| and the 
Kohn variational principle jl6|,^2]]. In this Letter, the 
wave functions of the scattering states are expanded in 
terms of the CHH basis |l4j] and the complex form of the 
Kohn variational principle is applied (2^]24|] . 

The wave function (WF) of a 4N state with total an- 
gular momentum J, parity n and total isospin T, T z (in 
the present case T — T z — 1) can be written as PI 
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The first term of Eq. ([j]) must be sufficiently flexible 
to guarantee a detailed description in the "internal re- 
gion," where all the particles are close to each other and 
the mutual interaction is large; goes to zero when the 
distance between the 3 He and the unbound proton i 
increases. This term in the WF is expanded in terms of 
CHH basis functions, following the procedure discussed 
in detail in Ref. p5| . 

The second term describes the asymptotic con- 
figuration of the system, for large Ti values, where the 
nuclear p- 3 He interaction is negligible. Let us introduce 
the surface functions 
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where \i is the spin function of the unbound nuclcon i 
and ^j^l is the 3 He bound state WF. This latter function 
is normalized to unity and is antisymmetric under the 
exchange of any pair of particles j, k and I. Hc has 
been determined as discussed in Ref. JTJ] by using the 
CHH expansion for a three-body system. 

The functions TZ L X \ri) of Eq. (||) are the ingoing (A = 
— ) and outgoing (A = +) radial solutions of the two-body 
Schroedinger equation without nuclear interaction p4| . 
The asymptotic WF is then written as 
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where the S-matrix elements '' 'S'^f, give the amplitude of 
the outgoing (L'S') component relative to the the ingo- 
ing (LS) wave. The elastic S-matrix, whose dimension- 
ality is 1 (2) for the J — (J > 0) states, should be 
unitary since there are no open reaction channels at the 
energies considered here. It follows that the eigenvalues 
of the S-matrix are written as exp(2i<5^ s ), where d"[ s is 
the eigenphase shift of the 2S+1 Lj wave. These quanti- 
ties are calculated by means of the complex form of the 
Kohn variational principle with a procedure similar to 
that one used in the N-d case |T4| , p4| |. 

Note that with the present method the unitarity of the 
S-matrix is assured only after the complete convergence 
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of the CHH expansion. Thus, for example, in the case 
of the J — waves, the calculated S-matrix will be of 
the form °<S££ = i] jn exp(2i(5° s ) where the "inelasticity" 
parameter t] JU may differ from 1 if the convergence is 
incomplete. For J > states, the inelasticity parameter 
can be defined as rj jn = ^Tr[S^ ■ S]/2. 

The expansion of the internal part ty^, is conveniently 
studied by grouping the functions of the basis in "chan- 
nels" (a given channel contains CHH states with the 
same angular-spin-isospin quantum numbers). The con- 
vergence of the L = waves (J n = + , 1 + ) at E cm . = 
was studied previously in Ref. [Q. At E c , m . = and 
at the energies considered here, a rather small number 
of channels is sufficient to provide a good convergence. 
This is due mainly to the Pauli principle which limits 
overlaps between the four nucleons. As a consequence, 
the internal part is rather small and does not require a 
large number of channels. 

In contrast, for the L — 1 waves ( J n = CP, 1~ and 2~) 
the convergence rate is slow and many channels have to 
be included. In these cases, the interaction between the p 
and 3 He clusters is very attractive (it has been speculated 
that some 4N resonant states exist) and the construction 
of the internal wave function is more delicate. 

An example of convergence for the 0~ and 2~ states 
is reported in Table | (in the J = 2 case the results 
arc relative to the L = S = 1 wave). The calculation 
has been performed using the AV18 potential at E c . m . — 
1.69 MeV for a few values of the number N c of channels 
included in the expansion of ^q. The values rj JU w 1 for 
N c = is accidental, in fact, the value of rj jn increases 
after including a few channels more. The convergence is 
reached only for N c 3> 10 when \rf m — 1| w 10~ 5 . 

At the energies considered here, the scattering in the 
L = 2 waves (J n = 1 + , 2 + and 3 + ) is very peripheral 
and the corresponding phase shifts are small. They can 
be calculated with good precision by considering only the 
asymptotic part in Eq. (|2|) . The contribution of the waves 
with L > 2 is very tiny and has been disregarded. 

The predicted analyzing powers are compared with the 
measurements in Fig. 0. The solid (dashed) curves corre- 
spond to the AV18 (AV18UR) interaction model. The 
main aspect to be seen in Fig. |] is that the calcula- 
tions significantly underpredict the analyzing power by 
approximately 30% at 1.20 MeV and 40% at IMMeV. 
Similar results have previously been seen and well doc- 
umented for N-d scattering. We also see in Fig. that 
the 3N interaction has almost no effect on A y . 

In Table || the theoretical phase shift parameters at 
E c . m . — 1-20 MeV are compared with those obtained 
in the energy-dependent phase shift analysis (PSA) of 
Rcf. p6[ which well reproduce the observables shown in 
Fig. 111. As can be seen, some phase shifts are well repro- 
duced by the theory. However, the 3 p2 and 3 Pi phase 
shifts are both sizeably underpredicted by the theory. 
For some partial waves the errors are too large to make a 



conclusive statement about the quality of the calculation. 

The A y observable is sensitive mainly to the P-wave 
phase shifts. At E c . m . — 1.20 MeV, for example, the 
A y maximum increases by approximately 42% when the 
3 p2 phase shift is changed by +10%. On the other 
hand, changing the 3 Po ( 3 Pi) phase by +10% decreases 
the A y maximum by 8 (6)%. This observable is in 
particular sensitive to the combination of phase shifts 
A = 0.5[ ( 5( 3 P 1 ) + ^( 3 P 2 )]-5( 3 P ). AtE c . m . = 1.20 MeV, 
the theoretical calculation predicts A = 3.9° (3.8°) with 
(without) the 3N force. The corresponding experimental 
result (from the PSA) is A = 6.9° ± 0.9°. 

It is interesting to note that the discrepancy between 
the theoretical and experimental A y is very much like 
that for the N-d case. There the main problem is that 
the splitting between the A P\/2 phase and the average 
of the 4 P3/2 and 4 Ps/2 phases is too small. For example, 
for p-d scattering at E c . m . — 2 MeV, the calculations give 
0.5[5( 4 P 3/2 ) + 5{ A P b/2 )] - 5( A P 1/2 ) = 1.87°, whereas the 
phase shift fits to A y require a splitting of 2.61° |J. It 
is plausible to suspect that the A y problems for N-d and 
p- 3 He both arise from the same deficiency in the nuclear 
Hamiltonian. Investigations in this direction are actively 
being pursued. 

At the energies studied here old differential cross sec- 
tion data exist jp7|. The theoretical estimates are com- 
pared with these data in Fig. 0. Overall the agreement 
is good, especially in the minimum region, but at large 
angles a small discrepancy is seen. This discrepancy can 
once again be traced to the undcrprcdiction of the calcu- 
lated 3 p2 and 3 Pi phase shifts. For this observable the in- 
clusion of the 3N force produces small but non-negligible 
effects. Finally, there are also some measurements of the 
3 He analyzing power in the energy range of interest p8[ . 
For this observable, the theoretical estimates are again 
below the data, but unfortunately the experimental un- 
certainties are too large to make a definitive statement. 

In summary, we have performed calculations of low- 
energy p — 3 He scattering observables, employing a real- 
istic interaction with a 3N force, and including the effects 
of P and D waves. A comparison with new proton an- 
alyzing power measurements demonstrates that the A y 
puzzle also exists in the 4N system. 

This work was supported in part by the National Sci- 
ence Foundation under grant number PHY-9722554. 
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TABLE II. S- and P-wave phase shift and mixing angle 
parameters (in degrees) at E c .m. = 1.2 MeV calculated with 
the AV18 and the AV18UR potentials. The values obtained 
with the energy dependent PSA of are also shown. 

wave 

3 Po 
3 Pi 
3 P 2 



FIG. 1. Measurements of the proton analyzing power A y as 
a function of the scattering angle at E c .m. = 1-20 MeV (panel 
a) and 1.69 MeV (panel b). The theoretical estimates ob- 
tained with the AV18 (solid curves) and the AV18UR (dashed 
curves) interaction models are also reported. 



FIG. 2. As in Fig. but for the differential cross section. 
The data points are from Ref. . 
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TABLE I. Eigenphase shifts S^s (in degrees) and inelas- 
ticity parameters T] jn for the p- 3 He scattering waves 3 Po and 
3 P 2 at E c . m . = 1-69 MeV, calculated with the AV18 potential 
and the complex Kohn variational method. N c is the number 
of channels included in the CHH expansion of the wave func- 
tions (the case N c = corresponds to including in the WF 
only the asymptotic part). 
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